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Notations 


k 

turbulent  kinetic  energy 

U,  u 

mean  and  fluctuating  velocity 

ke 

temperature  variance 

uv 

Reynolds  stress 

Prt 

turbulent  Prandlt  number 

V0 

normal  turbulent  heat  flux 

P 

pressure 

y+ 

dimensionless  coordinate 

q; 

heat  flux 

OCt 

turbulent  thermal  diffusivity 

R 

time  scales  ratio 

£ 

dissipation  rate  of  k 

Rt 

turbulent  Reynolds  number 

£e 

dissipation  rate  of  k@ 

t,  e 

mean  and  fluctuating  temperature 

X 

kinematic  time  scale 

t 

time 

xe 

thermal  time  scale 

Introduction 

The  very  high  temperature  level  reached  in  actual  turbine  components  requires  accurate 
simulation  tools  to  predict  the  heat  transfers.  The  turbine  flow  is  very  complex.  It  is  dominated  by  high 
pressure  gradients,  three-dimensional  viscous  effects,  and  high  heat  fluxes.  Moreover,  the  turbulence 
plays  a  major  role  in  the  heat  transfer  level.  Very  efficient  turbulence  models  have  been  developed  for 
the  prediction  of  dynamic  flows.  Modelisation  of  thermal  field  is  relatively  inferior  to  the  modelling  of 
dynamic  field.  The  current  practice  is  to  use  assumption  of  constant  turbulent  Prandtl  number.  The 
objective  of  this  work  is  to  develop  and  assess  thermal  flux  models  for  turbine  aerothermal  field 
prediction  [1], 


Physical  model 

The  physical  model  is  based  on  the  three-dimensionnal  compressible  Favre  averaged  Navier- 
Stokes  equations.  The  equations  are  written  in  a  conservative  form  for  cartesian  coordinates. 

Continuity  equation  : 

dp  !  dpt/,-  _  Q 

dt  dxj 

Momentum  equations  : 
dpU,  |  dpU,UJ  _  dPs  |  datJ  dpiijUj 

dt  d.Xj  dxt  dxj  dxj 

Energy  equation  : 

dpE  +  dpEUi  _dZjjU j  dq.  dpute 
dt  dxj  dxf  dxj  dxj 

where  Ui,  p  ,  P  and  p  are  respectively  the  velocity,  the  density,  the  pressure  and  the  viscosity. 
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These  equations  are  implemented  in  the  Navier-Stokes  code  CANARI,  developed  by  ONERA 
[2].  The  equations  are  discretised  in  space  with  a  central  finite  volume  scheme  associated  with  a  four 
step  Runge-Kutta  time  scheme.  Second  and  fourth  order  dissipative  terms  are  added  to  ensure 
numerical  stability.  An  implicit  smoothing  residual  could  be  used  to  accelerate  the  convergence. 


Turbulence  modelling 


To  close  the  momentum  equations,  the  Reynolds  stresses  m,m;  have  to  be  determined.  The 

classical  Boussinesq  hypothesis  is  postulated  which  links  the  unknown  Reynolds  stresses  to  the  known 
velocity  field  with  the  turbulent  eddy  viscosity  p,.  Based  on  a  dimensional  analysis,  the  modelisation 
of  p,  is  built  with  a  velocity  scale  and  a  time  scale.  The  turbulent  kinetic  energy  k  and  its  dissipation 
rate  £  are  used  to  express  these  two  scales  and  pt.  The  two  transport  equations  of  k  and  £  are  solved. 
Damping  functions  of  Hattori-Nagano-Tagawa  [3]  are  introduced  to  solved  the  equations  up  to  the 
wall. 


Thermal  flux  models 

The  turbulent  thermal  model  is  based  on  a  two-equation  model.  Similar  to  the  modelling  of  the 
Reynolds  stresses,  the  thermal  fluxes  are  linked  to  the  mean  temperature  gradient  by  a  scalar  :  the 
turbulent  diffusivity  at.  A  dimensional  analysis  is  earned  out  to  express  oct-  This  quantity  has  the  same 
dimension  as  pt.  The  most  representative  velocity  scale  is  that  of  the  turbulent  phenomenon,  so  the 
root  of  the  kinetic  energy  is  chosen.  The  time  scale  is  expressed  as  a  product  of  the  velocity  time  scale 
x=k /£  and  the  thermal  time  scale  Xe=ke/£e.  k@  is  the  temperature  variance  and  £e  its  dissipation  rate. 

=C-J\pk4r^ze 

This  kind  of  model  signifies  that  the  kinematic  and  the  thermal  phenomena  are  coupled  and  the 
k-£-ke-£e  equations  are  solved  simultaneously  [4],  [5],  [6].  The  ke-£e  equations  are  developed  for 
incompressible  flows  but  on  the  Morkovin  hypothesis  [7]  they  could  be  adapted  to  compressible  flows 
if  the  Mach  number  fluctuations  are  small  compared  to  the  mean  Mach  number  value. 

In  the  ke-equation,  the  only  term  to  be  modeled  is  the  turbulent  diffusion  and  a  gradient  type 
representation  (like  for  the  k-equation)  is  used.  Various  models  have  been  proposed  for  the  £e- 
equation.  The  model  we  implement  has  two  production  terms  and  two  destruction  terms.  For  the 
turbulent  diffusion  term,  a  gradient  approach  is  adopted. 

The  low-Reynolds  number  model  of  Hattori-Nagano-Tagawa  [3]  has  been  chosen  because  it 
solves  the  pseudo-dissipation  £0'  transport  equation,  which  allows  zero-value  at  the  wall.  The  use  of 

this  variable  implies  the  introduction  of  additional  terms  :  i^yjkg  /dxi )  and  Ee. 


The  modelling  equation  of  ke  is  : 
dpkd  dpUjkd  _  — -dT  d 

- 1 - —  — piljU - 1 - 

dt  dx,  dx,  dx 


a  +  - 


a, 


convection 


production 


LV 

diffusion 


dkH 


dx. 


P£a 


dXj 
dissipation 


The  modelling  equation  of  £e'  is  : 

dp£o  I  u  dp£»  -c  f  p£e  P  +  C  f  ^-P- 
dt  Ui  dx Pltp'~p^  6  p2fp2  k 


_a_ 

dx. 


{  a,  ^ 

d£e 

C  f  p£s 

^dlJd  1  - — 

LI  a*j 

dx, 

pe 2 

£ 

k 


convection 


production 


diffusion 


destruction 
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where  : 


fpi,fP2,fdi=  1,  o^,ae=  1,  P  =  ~puiuj^-L,  Pe  =-puid^L 


-dU, 
3x , 


37^ 

fir, 


/  =  [l-exp(-y+/30)]~ 


20 


1  + 


[-(R,/120)2. 


/*  = 


1-exp 


V  30 , 


;Ee=-a,(l-/J 


1  -  exp 


>’+^, 


1/3  \ 


30 


y-ii 


,  7.9 

1  + 
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P 

v  120  J 


32r 

3x,3x, 


=(l.9/£2-l)/0.9;  /£0  =(l-exp(-/P(K/3o| 


Channel  test  case 

This  test  case  is  a  wall-bounded  shear  flow  with  thermal  gradients.  The  results  of  the  calculation 
are  compared  to  the  Direct  Numerical  Simulations  of  Kasagi  [8].  A  constant  wall  heat  flux  is  applied, 
and  the  Reynolds  number  based  on  the  friction  velocity  and  the  half  height  of  the  channel  is  Re=150. 
The  buoyancy  effects  are  neglected.  The  grid  contains  200  points  in  the  axial  direction  and  60  points 
in  the  normal  to  the  wall  direction.  The  first  node  is  located  at  y+=0.05. 

The  mean  velocity  and  temperature  distributions  arc  presented  figures  1  and  2.  For  these  two 
quantities,  the  model  shows  a  good  behaviour.  In  particular  at  the  wall,  the  results  are  very 
satisfactory.  Figure  3,  the  calculated  Reynolds  stress  shows  very  good  agreement  with  the  DNS  data. 
The  same  remark  may  be  made  for  the  turbulent  heat  flux,  figure  4.  The  ratio  of  the  thermal  and 
dynamic  time  scales,  R,  is  presented  figure  5.  Its  evolution  is  qualitatively  reproduced  by  the  two- 
equation  model.  We  can  remark  that  R  is  not  constant  as  implied  by  a  constant  Prt  model. 


Turbine  blade  test  case 

A  turbine  stator  blade  is  now  tested.  The  experimental  facility  is  a  1  inear  transonic  cascade 
investigated  at  the  Von  Karman  Institute  [9],  [10]  .  The  Mach  number  is  0.151  at  the  inlet  and  0.85  at 
the  outlet.  The  stagnation  temperature  of  the  flow  equals  416. 8K.  The  wall  is  heated  at  Tw=296.25K. 
The  turbulence  intensity  at  the  inlet  is  4%. 

A  multi-domain  Fl-O-FI  grid  is  used.  The  upstream  FI  grid  contains  9x49x2  points.  A  O-type 
grid  is  defined  around  the  blade  to  ensured  a  good  description  of  the  wall  boundary  layers.  The  nodes 
distribution  is  301x65x2,  and  the  first  point  is  located  at  y+=  I .  The  downstream  FI-type  grid  contains 
1 89x49x2  nodes. 

The  heat  transfer  coefficient  FI  along  the  blade  is  available  and  compared  to  the  experimental 
data,  figure  6.  On  the  suction  side  (s/s0>0),  the  heat  transfer  coefficient  is  high  near  the  stagnation 
point.  Up  to  20%  of  the  curvilinear  abscissa,  FI  decreases,  due  to  the  effect  of  the  acceleration  on  the 
laminar  boundary  layer.  The  calculation  results  are  in  good  agreement  with  the  experimental  data.  At 
s/s0~0.25,  a  change  in  the  velocity  gradient  induced  the  beginning  of  the  laminar-turbulent  transition 
of  the  calculated  flow.  The  transition  produced  the  strong  increase  of  the  heat  transfer  coefficient.  In 
the  experiment,  the  laminar-turbulent  transition  appends  more  downstream  and  the  increase  of  FI  is 
delayed.  The  k-£  model  is  unable  to  simulate  precisely  the  transition  position.  When  the  transition  is 
completed,  FI  decreases  slowly  up  to  the  trailing  edge.  Along  this  region,  the  heat  transfer  level  is  well 
predicted  by  the  computation. 

On  the  pressure  side  (s/s0<0),  the  acceleration  of  the  laminar  boundary  layer  induced  the 
decrease  of  H.  At  s/s0=  -0.1,  a  small  separation  bubble  occurs  and  initiates  the  laminar-turbulent 
transition.  These  two  phenomena  produce  the  increase  of  the  heat  transfer  level.  Qualitatively,  the 
calculation  predicts  evolutions  in  agreement  with  the  experimental  ones.  But  the  transition  is  induced 
too  late  by  the  simulation.  Downstream,  the  heat  transfer  coefficient  increases  up  to  the  trailing  edge. 
This  increase  is  limited  by  the  effect  of  the  acceleration  in  this  region.  The  computation  seems  to  be 
unable  to  account  this  effect  and  overestimates  the  FI  level. 
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Concluding  remarks 

A  two-equation  model  for  the  turbulent  thermal  field  has  been  implemented,  in  addition  of  a  k-e 
model  for  the  dynamic  field.  These  four-equation  model  has  been  developed  in  a  Navier-Stokes  code. 
This  approach  is  more  universal  and  avoids  to  introduce  the  assumption  of  a  constant  turbulent  Prandtl 
number  which  is  not  valid  when  strong  temperature  gradients  exist. 

The  model  is  validated  for  a  channel  flow.  The  results  are  very  satisfactory.  The  heat  fluxes  are 
well  predicted.  A  turbine  cascade  flow  is  simulated.  The  physical  phenomena  are  reproduced  by  the 
model.  Although,  the  position  of  the  laminar-turbulent  transition  is  not  correctly  predicted  by  the  code 
and  induced  some  discrepancies  in  the  heat  transfer  level.  However,  the  capability  of  solving  this 
sophisticated  closure  model  in  a  full  Navier-Stokes  code,  for  complex  turbine  flows,  is  demonstrated. 
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